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Abstract 

We study a one-dimensional model of interacting conduction electrons 
with a two-fold degenerate band away from half filling. The interaction in- 
cludes an on-site Coulomb repulsion and Hund's rule coupling. We show that 
such one-dimensional system has a divergent Cooper pair susceptibility at 
T = 0, provided the Coulomb interaction U between electrons on the same 
orbital and the modulus of the Hund's exchange integral \J\ are larger than 
the interorbital Coulomb interaction. It is remarkable that the supercon- 
ductivity can be achieved for any sign of J. The opening of spectral gaps 
makes this state stable with respect to direct electron hopping between the 
orbitals. The scaling dimension of the superconducting order parameter is 
found to be between 1/4 (small U) and 1/2 (large U). Possible experimental 
realizations of such systems are discussed. 

PACS numbers: 74.65.+n, 75.10. Jm, 75.25.+Z 



1 Introduction 

The idea that a doped spin hquid state favours superconductivity (SC) belongs to 
Anderson |jl|]. This idea was introduced in the context of quasi-two-dimensional 
copper oxide superconductors which are believed to be adequately described by 
the t-J model. However, it has proved very difficult to reliably analyse this model 
on a square lattice. Recently it has been suggested]^, that SC might appear 
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in one dimensional doped spin liquids which are much easier to analyse [Q. In 
this letter we further develop these ideas and give a more detailed description of 
the SC spin liquid state suggested in these papers. 

We study the model of two Hubbard chains with the Hund's coupling between 
the chains: 

H = Huuh + XI {^Pl,rp2,r + 2 JSi,rS2,r) (1) 

Huub = -^Y. (c;t",r+lCa,a,r + H.C.) + U 4,aACr,a,^C^,a,l^r,a,l (2) 
r,a,a r,a 

where c^an Ca,a,r are electron creation and annihilation operators on the site r, 
a = ±1/2 and a = 1,2 are electron spin and chain (orbital) indices respectively, 
and 

Pa = Yl ct,aCa,a; ^ = ^Y. C^ta^^a/^Ca,/? (3) 

where cr" are the Pauli matrices. Electronic mechanisms give ferromagnetic ex- 
change (Hund's rule); an antifferromagnetic J can be generated by the Jahn- Teller 
effect (see the discussion later). Later we shall present the arguments that the 
model (1) may describe quasi-one- dimensional materials based on chains of Cgo 
molecules. 

The model (|1]) does not include a direct hopping between the orbitals because 
such terms are irrelevant when there are spectral gaps. We shall demonstrate that 
aX U » u such gaps exist for all excitations except for the symmetric charge 
density mode. Since a single electron hopping process inevitably involves modes 
with spectral gaps, the direct hopping becomes energetically costly and therefore 
irrelevant. It turns out that if the certain conditions on f/, u and J are met, the 
ground state of the model (1,2) has strongly enhanced SC fluctuations. Due to the 
presence of spectral gaps in a single pair of chains (further we shall sometimes call 
such pair a "ladder"), the only excitations which can tunnel between the ladders 
are Cooper pairs. As we shall demonstrate later, the scaling dimension of the SC 
order parameter is small which makes the pair susceptibility more singular than 
in BCS superconductors and thus increases the critical temperature of a three 
dimensional system made of weakly interacting ladders. 

Khveshchenko and Rice0 in their analysis of the model (1,2) have used 
bosonization in the weak coupling limit combined with the renormalization 
group approach (RG) (later this model has been also studied by Fujimoto and 
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Kawakami]^ with the same quahtative results). We suggest an alternative ap- 
proach which allows us to proceed directly to the strong coupling fixed point 
avoiding a cumbersome RG analysis of the weak coupling limit. The suggested 
reformulation of the Hamiltonian enables us to expand around the strong cou- 
pling point. The reason this representation is so useful is that it elucidates the 
inherent symmetries of the Hamiltonian. We demonstrate that in order to realize 
the mechanism of SC suggested in 0,0, one needs to have a hierarchy of cou- 
plings where the on-chain repulsion must be the strongest one. Together with 
the interchain exchange interaction it creates spectral gaps for all excitations ex- 
cept for the symmetric charge mode. The interchain Coulomb repulsion u has a 
potentially damaging effect on SC establishing a lower limit on the magnitude of 
\J\. It is intriguing that the conditions for SC can be met for both signs of J. 
The SC state is stable with respect to direct hopping if the hopping amplitude t 
is much smaller than a sum of the gaps. 



2 Bosonization and Introduction of Majorana 
Fermions. 

As usual, interactions between electrons in one-dimension cannot be treated per- 
turbatively. Hence we have to resort to non-perturbative methods, the most 
popular among which is bosonization. Since the Hubbard model is exactly solv- 
able for any U/D and doping, the bosonization approach requires only a smallness 
of the Hund's rule couplings u and J. The behaviour of the model at half 
filling has been already well studied In this paper we assume a fair amount 
of doping. 

We begin by bosonizing the charge and the spin density operators (^. Accord- 
ing to Frahm and Korepin(0) the continuous limit of these operators is given by 
(here we omit the orbital indices) 



p(x) = Jr + Jl + 2 cos(v^<I>s 



aee^^*^^^ exp(iy2^$c) + H.c. 
a^e^^'^^^ exp(iv^$c) + H.c. 



Six] 

= isin(\/2^$,)e-i^27r$ 



+ 

Jr + 3l + 



2\kpx 



aix) + H.c. 



= exp[iV27r(-<I), ± 6,)] 



(4) 
(5) 
(6) 
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Ho = I I dx 



where the corresponding bosonic fields are governed by the following Hamiltonian: 

(7) 

where the operators satisfy the standard commutation relations: 

[<^a{x),eb{y)] = ie{x-y)6ab (8) 

and Qc, ttu and are numbers which numerical values are known only in some 
limiting cases (we shall return to this issue later). In particular, a„ = for a free 
electron gas. The charge and spin currents Jr^l and Jr,l obey the corresponding 
Kac- Moody algebras (see, for example, [^). The constant and charge and 
spin velocities v^-, Vg depend on the interaction and doping such that 1/2 < 
Kc < I. The lowest limit is achieved for U = oo and Kc = 1 corresponds to 
the non-interacting case. The fact that the corresponding spin constant remains 
unrenormalized {Kg = 1) refiects the SU(2) symmetry of the individual Hubbard 
chain (0). 

Substituting expressions (^, ^, ^ into Eq.(|lD and keeping only non-oscillatory 
terms, we get the following expression for the interaction density: 

u 

Vint = V + —[Jri + JLi][Jji2 + Jl2] + '^J['iRl + Jli][Jr2 + Jl2] (9) 



V = 2cos(v^$c )['"kcP(cosA/47r$j" + cosA/47r$s ) 

+ J|asp(cosy47r$j" - cosViTT^^ + 2cos V47r6j)] (10) 



The 4fcir-scattering contributes the irrelevant operator cos 4-^7?$^ which we omit. 
Here we have introduced symmetric and antisymmetric combinations of the 
bosonic fields:$^^ = ($J^^ ± $^^^)/\/2 This transformation leaves the bosonic 
Hamiltonian (|^) invariant. 

The interaction of charge currents can be easily incorporated into the free 
Hamiltonian (^ by the change of K^s and velocities for the symmetric and 
antisymmetric charge modes: 



K. 



c 



1±' ^ 



'2TlVr 



-1/2 

Kc (11) 



We shall assume that K~ < 1, i.e. the interaction in the antisymmetric 
charge channel remains repulsive. Then the scaling dimension of the operator 
V dv = ^ + K~ < 2 is always smaller than the dimension of the product of spin 
currents, which is a marginal operator. Therefore we omit the latter term. In the 



4 



subsequent analysis we will take advantage of the fact that all bosonic exponents 
in the square brackets in Eq.(|TD|) have the scaling dimension 1 which means that 
they can be expressed as fermionic bilinears. The corresponding expressions in 
terms of Majorana (real) fermions have been derived in Ref. [^]: 



(27rao) ^cosV^G^ = i{RiLi - RoLq) , (27rao) ^cosV47r$s = i(-Ri-^^i + -Ro-^^o), 

(27rao)"' cos v^<l>+ = i(i?2i^2 + RsL^^U) 



where is the short distance cut-off and the fermions satisfy the following anti- 
commutation relations: 

{Ra{x), Rb{y)} = ^S{x- y)6ab] {La{x), Lb{y)} = ^6{x- y)6ab, {Ra{x), Lb{y)} = 0(13) 
Being recast in new terms the Hamiltonian (|^) becomes 



c 



3 

+ivs dx {Lad^La - RadccRa) 

a=0 •' 



(14) 



The interaction density acquires the following form: 



V = i cos V 47r$^ [mi (1 - 6afi)RaLa - ma-Ro-^o] 
mi = A'naQ[u\ac\ + J\as\ ), m2 = 47rao(3J|as| — u|ac| ) 



(15) 



If K~ < 1 (that is u/2'KV- < I - K^, which requires the Coulomb repulsion U 
to be stronger than u\) the interaction (|15D is relevant making the Majorana 
fermions and the field $~ massive. The modes a ^ and a = acquire different 
masses which we call Mj and Mq respectively with t standing for "triplet" . The 
triple degeneracy of three fermion modes reflects the fact that they realize the 
spin 1 representation of the SU(2) group 0. The mass of the bosonic field $~ 
is denoted M. The dimensional analysis gives Mt, Mq, M ~ Ji/(i~^c ). One can 
safely assume that the following averages do not vanish: 

(cosv^<l>,-) ^0; ^0 (16) 

In the limit K~ << 1 this can be rigorously proven by the perturbation expansion 
in K~ . Namely, one should rescale the field $~: — y/^^c ^^"^ expand the 
cosine term in Eg. ([T5|) :cos V AttKc^^ = 1 — 2TTKc{^^y + ... Then in the first 
approximation in Kc one obtains the quadratic effective Hamiltonian ([l7|). One 
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can check that corrections coming from higher terms are not singular and contain 
powers of the small parameter K^.. The appearance of the averages (|TBp brakes 
the discrete Z2 symmetry which is not forbidden in (1 + l)-dimensions. The 
corresponding order parameter is non-local in terms of the original fermionic 
fields! . 

3 Superconducting fluctuations. 

As follows from the previous discussion, the Hamiltonian for the field $^ remains 
unaffected such that this field is still a free bosonic field. The other fields become 
massive and their strong coupling limit can be qualitatively described by the 

effective Hamiltonian: 

3 

Hejf = i / dxVs [-Rad:,Ra + Lad^La + (Mt(l - 5a,o) " Mo5a,o)RaLa] 
■' a=0 

+\ I dx[v:{d.e-f + v;id,^:r + M{^:f] (i?) 

As we have mentioned above this Hamiltonian becomes exact in the limit 0. 

We have dehberately put a minus sign at Mq in Eq. ([T7|) . We shall see that the 
correlation length of the SC order parameter is infinite only provided Mf and Mq 
in Eq.(|l^) have the same sign. As follows from Eg . ([15|) , the masses are 

Mj ~ (M|a/ + J|a/)(cosv^$^); Mo ^ {3 J \a ^l"^ - u\a d"^) {cos V2^<^^) (18) 

Combining this criterion with < 1 we obtain the following criterion for the 
existence of an infinite correlation length for SC fluctuations in the model (1,2): 

K; < 1; (MlaclVla/ + J){3J - u\ac\'^ / \a,\^) > (19) 

The coefficients known only in the weak coupling limit U/D « 1 

where Oc = a^. It is also likely that \as\ » {ad close to half filling where charge 
fluctuations are suppressed. Thus we can resolve the inequalities (|19[) explicitly 
only in the limit of weak interactions where Kc ^ 1 — U/Anvc'- 

U >u] J > u/3or J < -u (20) 

Let ipR^L be right- and left moving components of the original fermions. Now 
we shall demonstrate that if the criterion (|T9|) is met, the susceptibility of the SC 
order parameter defined as 

~ expiv/i(<|.+ + e+){exp[iv/^(0; + $;)] ± exp[-iV^{Q; + $;)]} (21) 
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is singular. The choice of sign in the above expression depends on the sign of 
(cosv^47r$^) (broken Z2 symmetry). The order parameter includes exponents of 
the fields and 67 with scaling dimensions 1/4. In the Majorana approach it 
is most convenient to express these fields in terms of the order (a) and disorder 
(fi) parameter fields of the Ising model using the fact that the model of massive 
Majorana fermions describes a two dimensional Ising model off the critical point, 
where the fermionic mass is related to the deviation from T^. M {T — Tc). The 
corresponding operators are related as follows[||: 

exp(iv/7r<l>+) ~ 1X1H2 + io-i(T2; exp{-i^/7^Q^) ~ a^Ho - icxo/xs (22) 

Substituting these expressions into Eg. (pi]) we get 



A± = exp(iv/^0+) cos v^<I>;(/xi/i2 + iaia^) \ "^'^^ ^ A : e'^""^ : (23) 

where depending on the sign of {cos 

A = (cosy7r$7(cri(T20-3Ato))or(cosA/7r$7(/Xi/i2/i30-o)) (24) 

In the expression (^31) we have omitted the term proportional to sin \/7i^~ since 
the sinus has a zero average. As is well known, one of the operators a and fi has 
a nonzero average off the critical point. Which one does not vanish depends on 
the sign of (T — Tc), that on the sign of mass in the Majorana representation. In 
order to have A 7^ 0, the mass term of the zeroeth Majorana mode should have a 
sign opposite to the mass term of the other modes which gives as criterion ([19|) . 

It is easy to see that if M^Mt < the order parameter describes the charge 
density wave: 

0CDW{X) = Y.i'^RM'^LM ± ^R,2a^L,2a) (25) 



As follows from Eg . (|23D , at T = the correlation function of Cooper pairs has 



the following asymptotic behaviour at distances » M (M is the smallest gap): 
(A(.,.)A-(0,0)) ~ M^^^j^—^^ (26) 

with defined by Eg.(|TT]). Since Kc varies from 1 at f/ = to | at f/ = oo, the 
scaling dimension of the order parameter 1/4 < = < 1/2 and the Fourier 
transformation of the correlation function (pairing susceptibility) always diverges 
at g = 0. At finite temperatures there will be finite number of kinks interpo- 
lating between vacua with (cosV^$c) > and < 0. Thus the Z2 symmetry is 
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restored with the finite correlation length ^ ~ exp(Mfc/T), where ~ Mq/Kc 
is the kink's mass. This exponentially large correlation length will be completely 
overshadowed by the correlation length of the $+-field: ~ 1/^- 



4 Physical Realization 

One posible group of experimental systems where the described superconducting 
mechanism may work are materials based on chains of Cgo molecules. An isolated 
Ceo molecule in a crystal has a rich degeneracy and the Hund's interaction may be 
generated by local vibronic modes (the Dynamical Jahn- Teller effect; see Ref. |T0| 



for a review). In a structure where Ceo molecules are situated on well separated 
chains a two-fold local degeneracy may survive even in the presence of an in- 
termolecular hopping. Then the one- dimensional electronic band is split into a 
single degenerate band created by the s,Pz orbitals and a doubly degenerate band 
created by Px,Py orbitals. The measurements performed for RbCeo give many in- 
dications of a quasi-one-dimensional behaviour The direct band structure 
calculations done for A„Ceo-compounds with A = K, Rb show, however, that the 
resulting electronic bands have a three dimensional character [|12]. Nevertheless, 
even if these calculations are correct in the given case, one may expect that it 
would not be too difficult to produce truly one-dimensional CgQ-based materials. 



The authors of Ref. |12| remark that the interchain hopping is very sensitive to 
the distance between the chains and therefore the one-dimensional limit "might 
well be achieved by doping with substantially larger species" . These authors fur- 
ther suggest creating chemical compounds of Ceo with "spacer molecules" such 
as NH3. 

The model describing a one- dimensional chain of Jahn- Teller molecules with 
electrons belonging to two degenerate local orbitals hopping along the chain and 
interacting with local optical (vibronic) modes of the molecules was introduced 



by Manini et al. [|13| . The Hamiltonian of electron-phonon interaction is given by 

(27) 



r,n=x,y 



9Cr,a,a^ abCr,b,aUn[r ) duJr? 2 



where u^{r), u^lr) describe local vibronic modes. A single Ceo molecule has a 
degeneracy higher than 2, but it is assumed that such high degeneracy does not 
survive in a crystal. Neglecting all retardation effects related to the kinetic energy 
of vibrons and integrating over the vibronic modes we get the interaction (|l]) with 
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u = J = 2g^jk. As in the standard BCS theory, the retardation effects will give 
us a high energy cut-off of the order of the vibronic frequency = (fc/M)^/^. 



Santoro et al.'^^ have studied this interaction in the approximation where the 
Jahn- Teller bound states were considered first, the hopping was treated as a 
perturbation and the Coulomb interaction was ignored. Thus the dynamical 
nature of the vibronic modes appeared to be an essential feature of this approach. 
As a result they got the scaling dimention of the SC order parameter equal to 
1/4 in agreement with our calculations. The fact that the two approaches give 
the same infrared behaviour is remarkable, but is by no means unusual. This is a 
standard feature of one- dimensional systems where one almost always reproduces 
strong coupling results starting from weak coupling. This is due to the fact that 
the strong coupling appears eventually in the renormalization process. 



5 Conclusions 

A curious property of our results is that SC requires a strong Coulomb repulsion 
between electrons on the same orbital. It is this repulsion which leads to formation 
of the spin gap with subsequent freezing out some electronic degrees of freedom. 
There is also a charge gap corresponding to antiphase fluctuations of the orbital 
charge densities. In the situation where the exchange interaction is generated by 
the Jahn- Teller mechanism, this gap coincides with the spectral gap of vibrons 
at g = 2kp. 

The Coulomb interaction, however, plays a double role since the scaling di- 
mension of the order parameter increases when the Coulomb repulsion grows 
becoming equal to 1/2 at f/ = oo. The small scaling dimension of the SC order 
parameter increases the probability of having a relatively large transition tem- 
perature in a quasi-three-dimensional system built of weakly interacting chains. 
If if: is the interchain Josephson coupling, then the dimensional analysis gives for 
the temperature of SC transition the estimate 

Tc~Mo(t/Mor; r^=^^ (28) 

where Mq ~ (^i+i/C^^-^c ) is the spin singlet gap. In the conventional case d ^ \ 
and the transition temperature is exponentially small in Mq/I, but for the present 
model 1/4 < < 1/2 and 2/3 < < 1. 

We conclude this latter with a brief description of the physical properties ex- 
pected for a system described by model (1,2). First, the magnetic excitations are 
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spin triplets and spin singlets, as at half filling 0. As for a general 5=1 magnet, 
the triplet excitations have a gap (the Haldane gap) which we call M^. There is 
a singlet gap Mq as well; the corresponding excitations are spin and charge sin- 
glets and can be treated as RVB (resonance valence bond) excitations. Remaing 
gapless mode is the charge free boson which gives a dramatic enchancement to 
SC fluctuations or, if the criterion MfMo > is not met, to CDW. We expect a 
real SC transition in the quasi-three-dimensional systems. 
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